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ABSTRACT 

The general properties of heat exchangers are examined in 
this study. Mathematical models are developed for the mechanism 
of heat transfer that take place during 

1. Evaporation in two-phase flow 

2. Condensation 

3. Heating or cooling 

From these models, computer programs have been developed for the 
preliminary designing of the following heat exchanger tynes: 

1. Counter-flow cooler with no change of state 

2. Counter-flow heater with no change of state 

3. Counter-f low condensor 

^ . Counter-flow evaoorator with two-phase flow 

5. Cross-flow heater 

6. Cross-flow cooler 

These programs are designed to be used in the examination 
of the size requirments for heat exchangers used in the opera- 
tion of thermodynamic cycles. The cycles of interest in this 
study are ones proposed for the production of useful oower from 
sources not conventionally used at present. 
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I 



INTRODUCTION 



In recent years, many proposals have been made for methods 
to produce oower from sources not conventionally used. With the 
recent increase in the cost of fossil fuel energy sources, it 
can be expected that many such proposals will be forthcoming in 
the near future. 

One such proposal came from Professor Clarence Zener of 
Carnegie-Mellon University ( 1) . This proposed cycle would pro- 
duce power utilizing the’ temperature gradient in the tropical 
areas of the oceans. To do this, a working fluid with a high 
vapour pressure is heated by the water at the surface of the 
ocean. In the tropics, this is around 25°C. The working fluid, 
as a vapor, then drives a turbine to produce power. It is then 
condensed back into a liquid by the deep ocean water at a temper- 
ature of 5°C. 

Hilbert and James Anderson(2) in 1966 proposed a cycle simi- 
lar to that of Professor Zener. It relies on the heating of the 
ocean surface water by the sun as the means for powering the sys- 
tem. In 1966, the cost of such a power source was deemed unrealis- 
tic, therefore the project was discontinued. In the light of 
the current shortage of energy, this proposal might warrent re- 
examination. 

Another such proposal came from a husband and wife team. 

/Aden and Marjorie Meinel, optical physicists, proposed a differ- 
ent way of using the solar energy to produce unsable power(3). 
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Their proposed cycle would use solar energy collectors to heat 
a liquid metal. This liquid metal would heat the working fluid 
which drives a turbine to produce power. 

The use of geothermal energy, heat from the earth's in- 
terior, as a power source has also been suggested( 4 ) , This 
type of heat source would be used in a manner similar to the 
warm surface water in the Zener proposal. 

All of these proposed cycles require similar machinery. 

The heat source is used to heat or evaporate the working fluid. 
For this, a heat exchanger is needed. The working fluid then 
drives a turbine and is cooled or condensed in another heat ex- 
changer. A pump or compressor is used to raise the pressure of 
the fluid and the cycle begins again. 

The feasibility of these and other such proposals must be 
examined. One phase of this examination should be to obtain an 
estimation of the size of the equipment necessary for the opera- 
tion of the cycle. A major item, both in the performance of the 
cycle and the size of the equipment, is the heat exchangers that 
are required. 

It is proposed that this study be an examination of the re- 
quired heat exchangers. This will be carried out by the develop- 
ment of a computer program for preliminary design purposes. 
General properties of heat exchangers will be examined. Mathe- 
matical models will be developed for evaporation, condensation, 
heating and cooling heat exchanger design. 

The heat exchangers needed for the cycle proposed by 
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Professor Zener will be designed in this study. The author 
feels that this will be useful for two major reasons. First, 
it will serve as an example of how this study and the computer 
programs developed can be used. A careful examination of the 
Zener proposal can also be useful in determining if it is a 
feasible method of producing power. 
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II 



HEAT EXCHANGER PROPERTIES 



2.1 General Properties 

2.1.1 Introduction 

A heat exchanger is a device used to transfer heat bet- 
ween two fluids that are separated by a wall. It is the form 
of the separating wall and the flow pattern of the fluids that 
distinguish between types of exchangers. This study will examine 
shell and tube heat exchangers and tubular cross flow exchangers. 
Fins are attached to the walls of some heat exchangers to increase 
the heat transfer. The use and effectiveness of fins will also 
be examined. 

Cost evaluation is an important problem in the heat ex- 
changer design. A minimum cost, including intitial .cos t and the 
cost of pumping the fluids, is generally desired. The alterna- 
tives in cost are examined in this study. 

2.1.2 Shell and Tube Exchangers 

A shell and tube heat exchanger is made up of round tubes 
mounted in a cylindrical shell. The axis of the tubes are para- 
llel to the shell. Shell and tube exchangers are classified by 
the direction of flow of the fluids and by the method of reducing 
thermal stress between the tubes and the shell used (5 (page 18-3*0 )• 
Selection of the type of exchanger to be used depends on a variety 
of properties including cost, pressures and temperatures. 

Fixed-tube-sheet exchanger is the simplest to fabricate and 
has the least cost. In this type of heat exchanger the shell and 
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and tubes are rigidly held. Because of this method of connec- 
tion the thermal stresses are the highest of any type of shell 
anql tube exchanger. Another disadvantage of fixed-tube-sheet 
exchangers is the inability to clean the shell side of the ex- 
changer by mechanical means. 

To eliminate the high stresses of a fixed-tube-sheet ex- 
changer, an exchanger with an internal floating head can be used. 
With this type of exchanger the tubes can be removed for clean- 
ing. The cost of an internal floating head exchanger is high 
and there is much more of a leaking problem than with a fixed 
exchanger. Pot the same job, the two types of exchangers would 
have the same size requirements and the design method is the 
same. 

The most commonly used shell and tube exchanger is the 
fixed-tube-sheet type(6). For this study it will be examined. 

A factor that has a large effect on the design of a heat 
exchanger is the choice of fluid to be passed through the tubes 
and that to be passed on the shell side. A number of factors 
go into the making of this decision: ( 5(Page 18-35) ). 

1. Cleanabillty . The shell is difficult to clean and 
therefore should have the cleanest fluid. 

2. . Pressure. It is much easier and cheaper to made 

tubes to withstand high pressure than the shells. 

High pressure fluids should pass in the tube. 

3. Temperature. High temperature fluids cause high 
stress. This stress can be more easily allowed 
for in the tubes. High temperature fluids should 
pass through the tubes. 

4. Quantity. The best design can be obtained when the 
fluid with the smallest flow rate is placed in the shell. 
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A counter heat flow exchanger has the flow of the two fluids 
in opposite directions. Figure 2.1-1. In most exchangers the 
flow is not purely counter flow. In the entrance and exit re- 
gions of the exchanger, some cross flow exists. If the exchanger 
is more than twice as long as it is wide, this cross flow region 
is small and can be disgarded. 

Parallel flow exchangers have the flow of both fluids in 
the same directions. See Figure 2.1-1. As with the counter 
flow exchanger, the end regions have some cross flow which is 
disregarded. Parallel flow exchangers require more heat trans- 
fer surface area, for the same temperature rise, than counter 
flow. See Figure 2.1-2. Parallel flow also has the disadvan- 
tage, compared to counterf low, , that the fluid being heated can 
only have a temperature rise equal to fifty percent of the hot 
fluid temperature at inlet. Counter flow exchangers can, in 
the limiting case, have the fluid being heated reach the hot 
fluid inlet temperature ( 7 ) . For these reasons parallel flow 
exchangers will not be further examined in this study. 

2.1.3 Cross Flow Exchangers 

In a single pass cross flow heat exchanger, the two 
fluids move perpendicular to each other. Figure 2.1-1. 

Cross flow generally requires more heat transfer area than coun- 
ter flow heat exchangers for the same temperature rise. Figure 
2.1-1. Cross flow units are used in areas of special application 
as the exchanger can be made compact. These exchangers are most 
widely used in gas to liquid exchangers with finned surfaces (8). 

(More discussion of fins appears in section 2.1-. *0. 
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Heat Transfer Surface Area 




Cold Fluid Temperature Rise as a Percentage of 
Hot Fluid Temperature Loss 



Figure 2.1-2 

Comparison of Required Heat Transfer Area 
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The cross flow heat exchangers examined In this study will 
be tubular exchangers. The heating or cooling fluid is passed 
through a tube and remains unmixed. The working fluid flows 
across the tubes and is allowed to mix. 

2.1.4 Pinned Surfaces 

Pins are extensions of the tubes that increase the heat 
transfer surface. They are most effectively used when the heat 
transfer is between a liquid and a gas. A rough rule-of-thumb 
for the heat transfer surface area of the fin is that the sur- 
face area for the gas and liquid should be inversely proportion- 
al to the ratio of respective heat transfer coefficient ( 7(page 
182) ). This rule must be slightly modified because the effici- 
ency of a fin decreases with an increase in area. If, for ex- 
ample, an air-water heat exchanger was desired, the area of the 
air side should be 10 to 20 times as large as the water side. 
This is caused by the heat transfer coefficient of water being 
about 500 Btu/ hr ft °P to 20 Btu/ hr ft °P for air. To make 
up this need for different areas, fins would be used on the air 
side . 

The spacing of fins in a heat exchanger is generally de- 
termined by manufacturing limits. The height is limited by the 
required flow area. A spacing of from 6 to 15 fins per inch is 
the normal range (9). , 

In the heat exchangers examined in this study, fins are 
not used. Cost is an important factor. The fabrication costs 
of fins is much higher than for non-finned tubes. If the fin 
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is not either integrally or metallurgucally bonded to the tube, 
the fin efficiency may be adversely af fected(ll) . This requires 
the fins to be soldered or welded to the tubes. This is an ex- 
pensive fabrication process. Because fins decrease the flow 
area, they increase the required pumping power for the fluids, 
increasing the operational costs. 

The nature of the fluids used in this study also discourage 
the use of fins. In the proposed cycles, the heat transfer is 
from a fluid in one phase to a fluid of the same flow or to a 
two phase flow. Fin performance in both cases is less than ideal. 

If a detailed design of heat exchangers was going to be 
carried out, the use of fins should be included. As this study 
is a preliminary design to determine feasibility, examination 
of fins is not needed. 

2.1,5. Cost of Construction 

The estimation of cost of a heat exchanger is a diffi- 
cult task and is often only a very rough estimate. Many rea- 
sons for this exist. Companies that make heat exchangers and 
components generally consider cost of both materials and fab- 
rication highly proprietry. They are, therefore, quite reluc- 
tent to release this data. What available cost data there is, 
is a number of years old by the time it is published. With 
the current rate of inflation, data as little as six months old 
is inaccurate. 

Cost of a heat exchanger is sensitive to special require- 
ments. A small matter such as rigorous quality control can 
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easily double the cost of an exchanger ( 11) . The use of spe- 
cial materials, unusual size materials, or non-standard shapes 
can greatly change the cost of an exchanger. All of these fac- 
tors are difficult to include in a cost estimation. 

In this study a rough cost is examined. The values are 
from 1960(12), and are therefore quite inaccurate. They do serve 
to give a method of comparison between different shapes. If two 
exchangers are designed to do the same job, examining their com- 
paritive costs, the best can be selected, regardless of the ab- 
solute accuracy of the numbers themselves. 

Two cost factors are examined; the cost of the material 
making up the shell, and the tubes, and the cost of fabrication. 
The fabrication costs include header joint cost (welding), cost 
of drilling of the header sheet, cost of cutting tubes to size 
and positioning of the tubes. The cost does not include inspec- 
tion, sales costs, and shop overhead( 12) . The material costs 
are determined on a base size heat exchanger. Corrections are 
made to adjust the base size to the particular heat exchanger 
designed. The two correction factors used are the Tube Material 
Factor and the Tube Length Factor. The tube material correction 
takes into account the higher costs associated with the use of 
special materials. The correction for tube length decreases 
as the length increases. This accounts for the economy of scale. 
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Cost ■ CW + CTB 



where : 
and 



and 



CVJ = CVJV x N 

CWV = 3.0 if Dout is less than h in. 

CWV = 12/5 x (Dout - .5/12) + 3.0 if Dout is 
greater than h in. 

CTB = 15 (AREA - 100 ) x TLF x AREA x TMF 

5 -* 

tlf i.3 - (length - 8)^ 

10 
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NOMENCLATURE 



AREA 


Heat transfer surface area 


CTB 


Cost of materials to build heat exchanger 


cwv 


Cost of fabrication per tube 


CV. T 


Cost of fabrication of heat exchanger 


Dout 


Outer diameter of tubes 


h 


Heat transfer coefficient 


k 


Thermal conductivity of material 


L 


Length 


P 


Wetted parameter 


q 


Rate of heat transfer 


T . 


Temperature 


AT 


Change in temperature 


TLF 


Tube length factor 


TMF 


Tube material factor 


S ‘ 


Cross sectional area 



Subscripts : 



f 


With fin 


o 


Without fin 
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2 . 2 EVAPORATOR PROPERTIES 



2.2.1 Introduction 

Boiling heat transfer is the mode of heat transfer that occurs 
when a liquid changes phase to a vapor. There are two regimes 
of boiling, film boiling and nucleate boiling. The form that 
boiling will take in a given circumstance depends on the degree 
of superheat of the liquid. Figure 2.2-1 shows typical pool 
boiling data(l£!). Pool boiling occurs when the heating surface 
is submerged beneath the surface of the liquid. In this study 
the region of interest is where Tw - Tsat is between 5 and 15 
degrees, in the nucleate boiling region. 

Nucleation is the formation, growth and motion of bubbles. 

For nucleation to occur in a liquid, that liquid must be super- 
heated, its temperature raised above the saturation temperature. 

The nuclei is formed on a foreign object in the liauid. This 
object is usually a cavity on the heating wall but it can be 
suspended foreign matter with a non-wetted surface ( I ^ page 6). 

The bubble grows by heat conduction to the liauid-vapor inter- 
face. The size it will become is based on a balance between 
buoyant and surface tension forces. 

When boiling occurs in a confined region with force convect- 
ion present, two-phase flow results. It is this type of boiling 
heat transfer that is discussed in this study. A mathematical 
model for two-phase evaporation is developed. 
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Figure 2.2-1 
Pool Boiling Data 
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Radiation Coming into Play 



2.2.2. Two-Phase Evaporation 

Two-phase flow is a pattern that contains both liquid and 
vapor. It occurs when a flow is subject to forced convection 
in a confined space, a tube, annuli, parallel plates. The ob- 
jective of the designer of a two-phase evaporator is to determine 
the heat transfer and the pressure drop characteristics of this 
flow. From this data the size of the required heat exchanger 
can be found. 

The determination of the desired characteristics can be a 
complex task as these two properties are coupled thermodynami- 
callyOS, page *47). A change in the heat transfer changes the 
quality of the flow and the flow pattern. This changes the 
pressure drop. At the same time a pressure change alters the 
quality of the flow and therefore effects the heat transfer 
characteristics. 

The flow regimes change as heat is added&o, page 1*0. 

Figure 2.2-2. The flow enters as a subcooled liquid. Boiling 
will take place at the walls even though the liquid as a whole 
is below the saturation temperature. This forms a bubbly flow. 

At moderate velocities, greater than 3 feet per second, the 
bubbles will be fairly evenly distributed. When the liquid reaches 
the saturation point, bulk boiling will occur. As the quality in- 
creases the flow becomes annular with a thin layer of liquid on 
the walls and a vapor core. In this vapor core will be liquid 
droplets that will decrease in amount and size as the quality of 
the flow increases. At some point the heated surface will become 
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dry. This is called the burnout point. From here on the heat 
transfer mechanism changes and must be dealt with in a separate 
manner(5> page 13-35). 

With the complex flows and the coupling of the heat transfer 
and pressure drop properties, many of the solution techniques 
depend on a combination of analytical and experimental studies. 

A basis for any solution for a two-chase problem is the assump- 
tion that the flow is fully developed (i*?, page 48). 

2.2.3 Mathematical Model 

Two models can be used for two-phase flow(20, page 24). The 
'homogenous' model considers the two phase flow as a single phase 
with mean fluid properties. This model is also known as the 
'friction factor' model. The 'separated flow' model considers 
the phases to be segregated into two streams, one liquid and one 
of vapor. If the velocity of each of the streams is considered 
to be equal, the 'separated flow' model becomes the 'homogenous' 
model. As one would expect, the 'homogenous' model is valid for 
the flow of low quality with the separated flow more applicable 
for the higher quality annular flovr. For the mass flow expected 
in this study, G less than 5001b m /ft 2 sec , the separated flow 
model gives better agreement with experimental dataC!0> page 47). 
For this mathematical model the 'separated flow' model will be 
used. 

The following assumptions are made for the 'separated flow' 
model : Ql) 



(overleaf) 



1. The vapor velocity and liquid velocity are constant 
but are not necessarily equal. 

2. There is thermodynamic equilibrium between the 
phases . 

3. Empirical correlations can be used to relate the two 
phase friction multiplier <{> 2 and the void fraction a 
to the independent variables of the flow. 

It is this third assumption that allows one to separate the 
calculation of the heat transfer and the pressure drop. The coup- 
ling effect of the properties is taken into account by the empiri- 
cal correlations. 

To find the heat transfer in two-phase flow, the standard 
heat transfer equation q = A UAT is used. 

V/here: q = w 

Because there is a great deal of variation in the value of U 
over the range from a quality, x, of 0 for pure liquid to a quality 
of 1 for pure vapor, the heat exchanger is broken Into segments 
of equal Ax. 



is much less than h w and h^. 

To find h^, the local heat transfer of the heating fluid, 
the following relationships are used: 02 ) 



To find U: 



i=i+i+ 

U hft h w 



Pout - Din 
kmat 



The last term is generally disregarded as k ma t 



Pout - Din 



For Re < 2100 




1.86 k h 



x x 

D ' 5 



- 25 - 



For 2100 < Re < 10000 



hh = 



k h 

D 



(Re h » 67 - 125) 



For Re > 10000 

. 023 k h Pr h * 4 G h - 8 

D 1,2 p h » e 



These relationships are empirical equations based on results of 
many experiments. Other relationships to find h can be used. 

(10, page 192). 

To calculate h w , the local heat transfer coefficient for 
evaporation, with two-phase flow is a bit more difficult. Chen (23) 
postulated that there are two mechanisms which take part in the 
heat transfer process for boiling in two-phase flow. One is the 
macroconvedtive mechanism that is associated with bubble nuclea- 
tion and growth as takes place in pool boiling. Chen further 
postulated that these two mechanisms are additive to their contri- 
bution to the total heat transfer. Others have suggested this 
idea. RohsenowQ^) in 1952 and Bambill0.5) in 19^2 both concluded 
that the heat transfer effects of boiling and convective flow 
when occuring at the same time should be additive. 

For the contribution of the macroconvective mechanism, Chenft3) 
felt that the heat transfer coefficient h mac , should be calcu- 
lated in a manner similar to forced convection flow. 

(equation overleaf). 
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h mac * °-°23 p Re •' Pr-‘ 

This equation is modified to show the effect of the two-phases. 

>W ■ °- 02 3 Rei •* Pr x •* ^ F 

where: „ G w D (1 - x) 

— * * 



p . CPI m 
kl 



This leaves the only unknown as F. F is the Reynolds factor and 
is a function of X^. 



where : 

F = ( Re/Re x ) * 8 



and 





'1 - x' 


.9 


M 


,5 


V 




x J 




Sit 




.^VJ 



For experimental data a curve of F vs l/X^t is plotted. (Figure 2.2-3) 
An equation can be empirically derived. 



F = 0.15 



; i 

*tt 



+ 



2.85 1 

Xtt' 475 J 



The contribution of the microconvective mechanism is derived 
from the analysis of Foster and Zuber(2 6) on pool boiling. Foster 
and Zuber found: 

(equation overleaf) 
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Reynolds Factor 



- 28 - 



.33 



Nu b = 0.0015 (Re b ) ' 67 ?r 1 

= hb r b 
k l 



r b is the radius of the bubble and is derived from bubble buoyancy 
and surface tension forces. 



Where : 


rb = 


’AT kip 1 Cp 1 a' 


.5 , 


D l ' 




AP J 




[g APj 





Forster and Zuber in their analysis disregarded the fact that the 
degree of superheat is not constant across the liquid film on the 
wall material. For convective flow this is important as the 
temperature gradient across this boundary layer is made much 
steaper by the flow rate. 

Chen(2.3>) then rev/rote the Foster and Zuber equation to take 
this temperature gradient into account. 

h mlc . 00122 * r " cp£ pg Eg AT .» aP ,» s 



o s Vl - 2 ’ X' 2 'p, 



S is the supression factor and is plotted as a function of 
Re = Rei F ** 25 . 

The total heat transfer coefficient h w is the sum of h mac 
and h m i C . This correlation appears to give better comparison 
with experimental data than any other correlations 09, page 126). 

The heat transfer coefficient h w cannot be calculated in the 
above method for the region that is after the burnout point. 
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Figure 2.2-4 
Supression Factor 
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For this study it is assumed that burnout will occur at a quality 
of .95. For the low amount of superheat expected this appears to 
be a reasonable value. For the segiment of the heat exchanger 

from quality of .95 to quality of 1, the following should be 
used for h w : 

h w = 0.023 k v Pr v •*> 0 vr 
D 2 U V 8 • 

Once the value of U is determined, the required heat trans- 
fer area can be found for each segment , 

Area . -, Hf C Ax 
U AT 

Where : 

4T - T wall - T sat 

By assuming a geometry of the heat exchanger, the number of tubes 
and the diameter of tubes, the length of each segment can be 
calculated . 



Az = Area/irD n 

This completes the heat transfer calculations. 

To calculate the frictional pressure drop in two-phase flow, 
the Lockhart-Martinelli-Nelson modelUf, 27 ) will be used. Although 
this model was derived for adiabatic flow, it is valid for heated 
pipes (.5", page 1*J - 9). The heat exchanger is broken into the 
same segment as for the heat transfer calculations. The pressure 
drops are then summed for the total. 
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The relation of two-phase friction pressure drop to single- 



phase friction pressure drop is: 



Tap 


= [ap' 


l Az J 


TPF l Az J 



<f>v 



2 



V 



Where <J> V 2 can be found empirically as is plotted as a function 



of X tt . 


Figure 2.2-5. 






<J> V = 1 + 2.85 


(X tt °- 52 2 ) 


Then 


[ap] = f A? ] 

. Az Jtpf l Az Jl 


*1 2 



Where <f>i can be found empirically 

4*1 = 42- 
xtt 



f AP] = <p y 2 x [ap] 

(.AZJrppp (.AZjj^'Q 



Where : 



AP 



Az 



' LO 



♦09 G w us y v ,2 
8 P v D1 - 2 



if Re >2500 



And 



AP 



Az 



J LO 



- 128 G w by if Re < 2500 
8 P v ° 2 



A P 

By substitution , (—) can be determined. AP W is calculated by 

Az <ppF 

multiplying (AH) by the length of the segment Az. The sum 



‘ Az ' 



TPF 
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Figure 2.2-5 

Two-Phase Frictional Multiplier 
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of the AP W is the friction pressure drop in the evaporator. 

If this pressure drop is too large to be handled by the 
cycle in which the evaporator is to be placed, the geometry of 
the evaporator can be changed. Increasing the flow area should 
decrease the length and the pressure drop. 

This concludes the necessary mathematical model for two- 
phase flow evaporators. 
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NOMENCLATURE 



A 


Heat transfer surface area 


Cp 


Specific heat at constant pressure 


D 


Diameter 


Din 


Inter diameter of tubes 


Dout 


Outer diameter of tubes 


F 


Reynolds number factor 


G 


Mass flow velocity 


g 


Gravitational constant 


h 


Heat transfer coefficient 


h fg 


Latent heat of evaporation 


k 


Thermal conductivity 


n 


Number of tubes 


Nu 


Nusselt number 


P 


Pressure 


Pr 


Prandtl number 


AP 

Az 


Pressure drop per length 


AP 


Change in pressure 


q 


Rate of heat transfer 


r 


Radius 


Re 


Reynolds number 


S 


Supression factor 


T 


Temperature 


AT 


Change in temperature 
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NOMENCLATURE Ccont.) 



X 


Quality 


x tt 


Martinelli parameter 


Ax 


Change in quality 


w 


Mass flow rate 


z 


Length 


X 


Latent heat of vaporization 


u 


Viscosity 


p 


Density 


0 


Surface tension 


<p 


Two-phase frictional multilplier 



Subscripts 



b 


Bubble 


h 


Heating fluid 


1 


liquid 


lm 


Lograthnic mean 


mac 


Macroconvective 


mat 


Material of walls 


mic 


Micro convective 


sat 


Saturation 


TPF 


Two-phase flov; 


V 


Vapor 


w 


V/orking fluid 
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2.3 



CONDENSATION PROPERTIES 



2.3.1 Introduction 

Condensation is the removal of heat from a vapor to convert 
it into a liquid. This is generally done by using a surface with 
its temperature T w , lower than the saturation temperature, T sat , 
of the vapor. The vapor will condense into a liquid on the colder 
surface. The exact nature of this condensation mechanism is not 
fully established. It is felt that the condensation begins in 
small cavities in the surface much as bubbles are formed in boil- 
ing(S", page 12-2). 

There are two forms of condensation - dropwise and film. 
Dropwise condensation occurs on a non-wetting surface. Droplets 
form and run down the surface, increasing in size by further 
condensation and coalescence. New droplets are formed to take 
their place. The heat transfer for dropwise condensation is 
from 5 to 50 times that for film condensation (31) . It is diffi- 
cult to maintain dropwise condensation for more than a short 
period of time without treating the heat transfer surface. 

Film condensation occurs when the liquid being formed wets 
the surface and establishes a stable film. This type of conden- 
sation is the prevalent in commercial condensers. As the heat 
transfer rate is lower with film condensation, condens rs are 
designed for this form of condensation. If the condensing should 
be dropwise for part of the operation of the condenser, it will 
perform better than design. In this study, film condensation 
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will be examined 



2.3.2 Film Condensation 

The first attempt to analyse film condensation was carried 
out by Nusselt in 1916. When the stable film is formed, vapor 
condenses on the liquid. The condensation takes place at the 
liquid-vapor interface because the heat is being transferred 
through the liquid. At this interface there is a continual 
interchange of molecules being condensed and molecules evapora- 
ting, with the net flow condensing. 

In his analysis Nusselt made the following assumptions: 

(10, page 31*0. 

1. The flow of condensate in the film is laminar. 

2. The fluid properties are constant. 

3. The condensate temperature distribution is linear. 

4. There are no changes of momentum through the film. 

5. The vapor is stationary and there is no shear force 
at the liquid-vapor interface. 

6. Heat transfer is by conduction only. 

These assumptions are important in the later development of the 
mathematical model for condensation. The only modification to 
these assumptions is in the temperature distribution. The actual 
temperature distribution is slightly curved(lO, page 239). 

2.3.3 Mathematical Model 

In the design of a condenser, the desire is to determine the 
required heat transfer area. The rate of heat transfer, q, is 
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found by: 



q = U AAT lm 

Where : 1_ a 1 + _1_ + Pout - Din 

U h c ^mat 

As was done in boiling heat transfer, K mat /Dout - Din can be 
disregarded. 



q = w h 



fg 



and 



ATi = ~ ^sat ~ ^c2) 



In 



T . - T . 
sat cl 



^sat “ ^c2 



V/ here AT^ m is the ' logarithic* mean temperature. 

To find the required area, all that needs to be calculated 
is the overall heat transfer coefficient, U. The heat transfer 
coefficient for the cooling fluid, h c , is found by the following 
relations (20) . 

if Re c < 2100 

h c = l,86(Re c •= Pr c ' s 2l) 

if 2100 < Re c < 10000 

h c = .116 — (Re * 67 - 125) Pr c * 67 
D c 

(cont. overleaf) 
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if Re c > 1000 



h 



c 



0,023 kcPrc 4 Gc 8 

y c 8 D 2 



To derive the coefficient of heat transfer for condensation, 
h w , one begins with the Nusselt analysis. For condensing on 
a verticle plate, 



Where h£ g = h fg + .68 Cp AT. 

The deriviation of this equation is examined in detail in References 
S', 31 and/C'. 

This coefficient must be modified for condensation on a hori- 
zontal tube. The relation above is valid with the modification 
that g sin is used to replace g. The average value of h for the 
tube is found by integrating h w from 0 to 



For a number of tubes in a vertical row, all of the conden- 
sate dropping from any tube is assumed to fall on the next lower 
tube. This will increase the film thickness of lower tubes and 
require a modification in h w . 



4 

h w = 0.943 



gPl(Pl ~ Py) 3 hfg 



Lp x AT 




\| Doutui AT 




\j n^ Dout AT 
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n. is the number of tubes in a vertical row. Observations by 
b 

Grant and OswentG'i) show that the condensate seldom falls from 
the upper tube to lower tubes as a continuous sheet. This would 
tend to lessen the importance' of having the n^ term. 

Two other effects can change the heat transfer coefficient. 
These are the presence of noncondensable gas and changes in the 
surface geometry. 

The presence of even a small quantity of non-condensable 
gas in a condensing vapor has a profound effect on the heat 
transfer coefficient. The non-condensable gas tends to form at 
the liquid-vapor interface. This increases the partical pressure 
of the gas at the interface and increases the resistance to heat 
transfer. The details of the methods to calculate this effect 
are given in Reference 31 . In this study, the effects will be 
disregarded. In most commercial condensers, ejectors are used 
to remove these gases and alleviate the problem. 

Changes in the surface geometry of the tubes can greatly 
increase the heat transfer coefficient with film condensation. 

One of the most promising of these geometry changes is fluting 
the tubes. This geometry was first examined by Gregorig(3i) . 

In experiments Gregorig obtained an h for condensing of 
8000 Btu/hr. ft 2 °F, a factor of four above what is obtainable 
with round tubes. Because of the complex shape of these tubes, 
it is difficult to analytically evaluate the heat transfer co- 
efficient. Instead, experimental data must be used. 

Upon finding h c and h w , the overall heat transfer coefficient 
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can be found and the required heat transfer surface area 
determined. A heat exchanger geometry, tube diameter and number 
of tubes is chosen by the designer. The condenser length can 
be found from this. 

1 = Area/im Dout 

The frictional pressure drop of the cooling fluid can be 
calculated. 



if Re > 2500 
c 

_ 0.046 *1 1 Oc* 

Re 0 - 2 D 2g p c 

Din G c 

V/here Re_ = 

• u 
M c 

if Re < 2500 
c 

P = _§i!_ 4 1 °c 2 

Re c Din 2g p c 

If this pressure drop is too large for the desired application, 
the condenser geometry, tube diameter and number of tubes, can 
be changed to decrease the pressure loss. 

This concludes the necessary mathematical model for con- 
densing on the outside of tubes. Condensing on the inside of tubes 
is not often carried out as it involves two-phase flow and is diffi- 
cult to analyse. 



NOMENCLATURE 



A 

Cp 

Din 

Dout 

g 

G 

h 

h fg 

k 

L 

n 

nb 
Pr 
AP • 

q 

Re 

T 

AT 

U 

w 

P 

P 

Subscripts 

c 



Heat transfer surface area 

Specific heat at constant pressure 

Inner diameter of tubes 

Outer diameter of tubes 

Gravitational constant 

Mass flow rate w/Area 

Heat transfer coefficient 

Latent heat of condensation 

Thermal conductivity 

Length 

Number of tubes 

Number of tubes in vertical bank 

Prandtl number 

Change in pressure 

Rate of heat transfer 

Reynolds number 

Temperature 

Change in temperature 

Overall heat transfer coefficient 

Plow rate ' 

Absolute viscosity 
Density 



Cooling fluid 
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NOMENCLATURE ( cont . ) 



1 


Liquid 


lm 


Logarithmic mean 


mat 


Material of the walls 


sat 


Saturation 


V 


Vapor 


w 


Working fluid. 



2.4 Cooling and Heating Properties 



2.4.1 Introduction 

In this study, the evaluation of the properties of cooling 
a working fluid and heating a working fluid in a heat exchanger 
are carried out at the same time. Thermodynamically the mechan- 
ism of heating and of cooling is the same. Any differences that 
do exist in the relationships will be discussed in the develop- 
ment of the mathematical model. 

As was noted in section 2.1, counter— flow heat exchangers 
and cross-flow heat exchangers will be examined in this study. 

The heat transfer mechanisms of both exchanger types are similar 
but the approach to the designing of each is different. 

The job to be done by these heat exchangers will be to trans- 
fer a specific amount of heat from one fluid to another. If the 
exchanger is a heater, the working fluid will be the cooler fluid. 
For a cooler, the working fluid must be the hotter fluid. The 
mass flow rate, w, of both fluids in each type of exchanger will 
be specified as will be the inlet and outlet temperatures. This 
will be true regardless of the flow pattern of the fluids. 

2.4.2 Counter-Flow Heaters and Coolers 

The mathematical model for the sizing of counter-flow 
heaters and coolers is based on the solving of similtaneous equa- 
tions . This is accomplished by setting groups of known constants 
equal to a large constant, Kl, K2, K3, etc. The heat transfer 
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area, number of tubes, and the exchanger length will be deter- 
mined in terms of these constants. This method is described in 



detail in Reference 7. 

To make this model valid for both heating and cooling, the 
hotter fluid in the exchanger is designated by a subscript 1 and 
the cooler fluid by a subscript 2. In a heater the working fluid 
is colder and is designated by the subscript 2, while for a cooler 
the inverse is true. The mathematical relationships are otherwise 
identical. 

A limiting constraint on the design of a counter-flow ex- 
changer is the pumping power required to overcome the friction 
pressure loss of the fluid flowing in both the tubes and the 
shell. In this design method, this constraint is set by estab- 
lishing an allowable pressure drop, AP. 



If Re! > 2500 




G 1,8 L 



(35) 



g P 1 J 



Din 1,2 




If Re, < 2500 



s 128 y^' 2 0 j 1 ^ L 



g p ]_Re r 8 j Din 1,2 



= K1 G-l 1 - 8 L 



Din 



1.2 



If Re 2 > 2500 

p a f °-092 Mg ' 2 ) G2 1,8 L 

g p 2 J Deq 1 * 2 

n 1.8 r 

P 2 = K2 _2 " 

Deq 1,2 



If Re 2 < 2500 







P 0 = 


fl28 y 2 2 ] 


G2 1 ' 8 L 


= 


K2 


G 2 1 * 8 L 






d 


g P 2 Re/j 


Deq 1 - 2 






Deq lwZ 


Where : 


Deq 


= (1.271 (s/Dout) 2 - 


1) Dout 




(7, 


page 324) 


This is 


for square pitch exchangers with s 


the 


spacing between 


tube centers. 














Where : 


Re ^ 


= Din (w^/ 


(ir Din2 ) 

k 










and 


Re 2 


= Deq (w 2 / 


Or Dec1 - 2 ) 
4 










Next the 


heat 


transfer coefficient, h, is 


calculated. (30, page 




If 


Re^ > 10000 














h 0.023 ki Pri 4 Gi' 8 


= 


K3 


Gl* 8 








Pi 8 


Din 2 






Din' 2 




If 


2500 < Re x 


< 10000 














hi - ' 


,116 kl Prr 33 Gi' 8 


= 


K3 


Gr 8 






Re^ 13 


Din’ 2 






Din 2 




If 


Rex < 2500 
















, 1.86 ki Prr 5 G !•» 

h l “ 


= 


K3 


Gr 8 






V 1 Re l 


Din- 2 






Din- 2 
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If Re 2 > 10000 

0.023 k 2 Pr 2 ’ 4 



h 2 = 



.8 



G 2' 8 



Deq' 



,•2 



.8 



= K 4 HL 



Deq" 



If 2500 < Re 2 < 10000 
h 2 



.116 k 2 Pr 2 * 33 



Re 13 
2 



G 2’ 8 = G 2 



Dea 



.8 



Deq 



If Re 2 < 2500 

h = 1.86 k 2 Prg- 



y 2 8 Re 2 3 



Dea : 



, Go - 8 
2 



Then: — ^ (both known) 

A Po K2 



f G l 



G 2 



1.8 



Deq : 



f Dea 'i K2 

[DinJ 



By definition _ w^/Ax^ 

g 2 w 2 / Ax 2 

Where Ax is the cross sectional area of the flow. 

Ax^ = tt/1( Din 2 n 

Ax 2 = tt/1( Dea 2 n = n/Jj Deq Dout n 



By substitution 



£i = 


r w ^ Dout' 


.4 


K2 K5 


G 2 


Wj bin ^ 




K1 



.33 _ 



= K6 



G 1 = W 1 Deq Dout 
G 2 Wj Din 2 



Then 



From solving for AP1/AP2 and G3/G2 



_ K3 


'2i\ 6 


[Deal 


.2 






G 2 j 


[DinJ 






_ K3 j 


\K2 K5 


.33 


DoutP 


f— i 


K4 1 


i K1 J 




LDin J 


IW2J 



K7 



The overall heat transfer coefficient can be calculated 

by: 

— = — + + Pout - Din 

u h l h 2 k mat 

The last term can be disregarded as it is insignificant compared 
to the other terms. This equation does not include the effect 
on U of scale deposits. For a preliminary design it can be dis- 
regarded. 

1 = _ 1 _ + _ 1 _ 

U h]_ h 2 

1 _ 1 Din' 2 J_ K6 Deo’ 2 K8 Din’ 2 

U K3 G^ 0 G^ 8 Gj 8 

Where K8 = -j-v K6 Dln ' 2 

K3 Kh Deo' 2 



The heat transfer rate 



q = A U A T im 
Where A = n Din L 

it • , 

q = n mDin U L A T^ m = ^ Din 2 n G^ Cp^ AT^ 
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Where AT-^ is the difference between the inlet and outlet tempera- 
ture of the hotter fluid and AT^ m is the logarithmetic mean 
temperature. 

Then 



L = £P Ti Gj Din 5 K q 



Din 



^A T 



lm 



U 



U 



L = K9 K8 D-l 1 - 2 G 1 



1.2 



Then by substitution of 



gt 



.2 



Din 



jj K9 K9 



Into 



K1 L G, - 8 
APn = L_ 



Din 



1 2 



Gt - 



L K1 K8 K9J 



.5 



= K10 



By definition 



G]_ = w^/Axi 



Therefore 



n = 






G 1 77 D l 2 



The value of Din and Dout are selected so L can be calculated, 



L = K9 K8 Din 1 ' 2 K10 



and the heat transfer area can be found, 

A = Din n L 
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This concludes the necessary calculations for the heat exchanger 
design. This design method assumes that: 

(1) The specific heat of the fluid, Cp is constant. 

(2) The overall heat transfer coefficient is constant. 

If these factors are not constant, the design aoproach is 

modified. The above procedures are carried out using average 
values of the fluid properties. A value of n and is determined. 
The exchanger is then divided into segments such that Cp and U 
can be assumed to be constant over the segiment. The design method 
is then repeated for each segiment with n and Gi kept constant. 

The length is determined for each segiment. These lengths can then 
be added to get the overall length of the heat exchanger. 

2.4.3 Cross-Flow Heaters and Coolers 

In addition to the information about the job to be done 
by the heat exchanger noted in section 2.4.1, one must know the 
condition of the flow through the exchanger to size a cross-flow 
heater or cooler. (3, page 18 - 8l). As stated in section 2.1.3, 
for this study it is assumed that the working fluid is flowing 
across the tubes and is mixed. This keeps the working fluid at 
a uniform temperature as it leaves the exchanger. 

In the development of the mathematical model, the equation 
q = U AAT m is again used with a slight modification - 

q = U A FAT lmc : 

(cont. overleaf). 
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Where 



(T - T ) - (T, - T „ ) 

Am — n2 c i. nJ. c2 

a X -i 



'lrnc 



In T h2 - T cl 



T hl " T c2 



A -lmc ^- s the l°garithmetic mean temperature for counterflow arrange- 
ment. This is not the integrated mean temperature AT m . F is then 
defined as: 

F = AT m /AT lmc (30) 

and is a function of e c and Z. 



Where: e c is called the effectiveness of the heat exchanger, 



and 



e c - T c2 - T cl /T hi - T ol 



Z = T hl ~ T h2 
T ol - T c2 



(30) 



F can be found on tables as a function of Z and e Q . F also 
varies with the flow pattern, if one or both of the fluids are 
mixed. These tables can be found in References 5, 10 and 30. 

The value of a can be found from a = w Cp^ ATj. All of 
these values are known. To find the heat transfer surface area: 



A = q/U PAT lmc 

U can be calculated by the relations for h^ and hg examined 
in the section 2.^.2 of this study. 

By selecting an exchanger geometry, number of tubes and tube 
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diameter, the length of the required exchanger can be determined. 

L = Area/uDin n 

The frictional pressure loss is then calculated. The relation- 
ships are examined in detail in section 2. *1.2 and need not be re- 
derived. If the pressure loss is too great, the geometry of the 
exchanger can be modified and the length recalculated. 

This design method assumes that U is constant throughout 
the length of the exchanger. If this is not a valid assumption, 
the exchanger can be broken into segiments. The above design 
procedure is then used for each seglment with the sum of each 
segiment length the total length. 
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NOMENCLATURE 



A 

Ax 

Cp 

Din 

Dout 

Deq 

F 

6 

G 

h 

k 

K 

L 

n 

AP 

Pr 

q 

Re 

s 

T c 

T h 

ATim 

U 



Heat transfer surface area 

Cross sectional flow area 

Specific heat at constant pressure 

Inside dianeter of tubes 

Outside diameter of tubes 

Equivalent diameter for shell side 

Mean temperature difference factor 

Gravitational constant 

Mass velocity 

Heat transfer coefficient 

Thermal conductivity 

Constant 

Length 

Number of tubes 
Frictional nressure drop 
Prandtl number 
Heat transfer rate 
Reynolds number 
Tube spacing 

Temperature of cooler fluid 
Temperature of hotter fluid 
Log mean temperature difference 
Overall heat transfer coefficient 



NOMENCLATURE (cont.) 



w Mass flow rate 

Z Temperature coefficient 

e c Effectiveness of heat exchanger 

y Viscosity of fluid 

p Density of fluid 



Subscripts 

1 Hotter fluid 

2 Cooler fluid 

mat Tube material 



- 55 - 



2 . 5 Material Selection 



2.5.1. Introduction 

In the selection of material for the tubes and other 
heat transfer surfaces in the heat exchanger, it is desirable 
to have the heat transfer coefficient as high as possible. 

An increase in the heat transfer at the surface of the material 
produces a corresponding increase in U, the overall heat trans- 
fer coefficient. This will decrease the reauired heat transfer 
surface area as the area proportional to the inverse of U. 

This augmentation of the heat transfer can be accomplished 
in many ways. The surface of the tubes can be roughened by a 
number of methods. The heat transfer surface can be left the 
same but the flow disturbed. This increases the heat transfer 
by displacement of the fluid or by setting up a vortex flow. 

Other methods of augmentation include vibration of the heat 
transfer surface, fluid vibration, and the use of electrostatic 
fields. Of these methods of increasing the heat transfer, the 
best appears to be the roughening of the surface material. The 
vibration and electrostatic methods require outside power to 
operate and are therefore less attractive. ( 5(page 10-2) ). 

2.5.2. Methods of Enhancement 

The first method for enhancing the heat transfer to be 
seriously considered, was increasing the surface roughness of 
the transfer material. The most successful for single-phase 
flow appears to be a method devised by Dipprey and Sabersky ( 1*J ) , 
and one developed by Kemeny and Cyphers (15). Dipprey and Sabersky 
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produced the roughness by electronically depositing copper on 
a rod coated with sand grains. The rod was dissolved leaving 
a tube of copper with close-packed sand grains on the inner 
surface. This method, for high Prandtl number flow, increases 
the h, heat transfer coefficient, by a factor of two. Figure 
2.5-1. The Kemeny and Cypher method uses a helical groove in 
the material and a helical protrusion from the material. The 
grooved surface showed little improvement in the heat transfer. 
The helical protrusion did increase the h by up to a factor of 
two. Figure 2.5-2. As can be seen from Figures 2. 5-1,2, the 
effects of the enhanced surfaces decrease with the increase in 
the Reynolds number. These methods also show no advantage in 
two-phase flow. 

In two-phase flow, the heat transfer is increased by pro- 
moting nucleate boiling at as low a heat flux as possible. 

Young and Hummel (16) placed spots of teflon on the heated sur- 
face. This promoted nucleation as shown in Figure 2.5-3. This 
lowered the required level of superheat required to produce 
boiling. With an increase of Reynolds number, the effect of the 
teflon spots decreases. 

The methods of increasing the heat transfer can also in- 
clude the use of a smooth surface with a higher k. All of these 
methods of increased heat transfer would be useful in the design 
of the heat exchangers for the thermal cycles . .What is needed 
is a method of comparison of these increases heat transfer sur- 
faces to unenhanced surfaces. 
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Figure 2.5-1 

Dippry and Sabersky Results 




Figure 2.5-2 

Kenery and Cypher Results 
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2.5.3 Methods of Comparison 



A method of comparison has been developed by Professor 
W.M. Rohsenow and others working in the Heat Transfer Lab at 
M.I.T. (17). The data on the surfaces is presented based on the 
same heat transfer area and tube diameter of equivalent diameter. 



f Re 3 is plotted agains-t j Re. Figure 2.5-3. 



Where: f Re 3 = P/v 2gp 2 De 3 a/A/V y 3 



and 



Ah Pr ^ 3 Pea 
V A/V Cpy 



The following are kept constant: 
g> p, Deq, A/V, y, Cp, Pr 



Therefore : 

f Re 3 ~ P/V 



and 



J 



Re 



Ah Qr n/ATi m 
V V 



where P = WAP/p 

For the same mass flow rate and material properties: 



A h/V ~ 



Ah 

WCp 

V 



NTU 

V 



For any flow arrangement there is one NTU vers e curve (18). 
Figure 2.5-4. 



Where : 



ATrise in cold fluid 
Thot in - Tcold in 
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f Re^ 

Pigure 2.5-3 



£ 



NTU 

Figure 2.5-4 
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The two materials can now be compared at two points. See 
Figure 2.5-3. The mass flow rate, w, the Deq, A/V, T hofc j_ n and 
T C old in remain constant throughout this comparison. 

At point A, tha values of w,' V, and Re are the same as at 
point 0. This gives the same size and shape heat exchangers with 
different surface material. The amount of heat transfer and re- 
quired pumping power will be different. These can be plotted as 
shown in Figure 2.5-5. 

Point B has the same mass flow rate, pumping power required 
and volume as point 0. A heat exchanger operating at B will have 
the same volume as one at 0, but the shape changes and there is 
more heat transfer. This comparison is made in Figure 2.5-6. 

If the comparison of the materials is to be made while doing 
the same job, point C and point 0 are examined. At these points 
the same heat transfer (NTU), mass flow rate and pumping power 
occur. The difference is that less volume is needed with materi- 
al B (point C). Figure 2.5-7. 

Since A/V has been held constant, the area required at point 
C is also less. This can be an important comparison if the heat 
exchanger needs to be as small as possible. A cost comparison 
can easily be done. As the pumping power is the same, the only 
variable cost is the material costs. 

Cost n / Souare foot x A n 

O D in, 

= Cr 

Cost./ Scuare foot x a 
A * A 
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Figure 2.5-5 

Comparison of Exchanger at Point 0 and 
at Point A 




Figure 2.5-6 

Comparison of Exchanger at Point 0 and 
at Point A 
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If Cr is greater than 1, A will be a cheaper material 
to use. If, however, Cr is less than 1, B is the cheaper 
material. 

A fourth comparison can be made at point D. Here the heat 
transfer, mass flow rate and the volume are kept constant with 
point 0. The shapes differ and the pumping power is lower. 
Figure 2.5-8. If the cost of the required pumping power was 
too high with material A, material B could be used to lower it. 
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Figure 2.5-7 

Comparison of Exchanger at Point 0 and 
at Point C 
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Figure 2.5-8 

Comparison of Exchanger at Point 0 and 
at Point D 
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NOMENCLATURE 



A 


Heat transfer surface area 


C P 


Specific heat at constant pressure 


Cr 


Cost ratio 


Deq 


Equivalent diameter 


f 


Friction factor 


g 


Gravitational acceleration 


h 


Local heat transfer coefficient 


k 


Thermal conductivity 


NTU 


Number of exchanger heat transfer units 


P 


Pumping power 


AP 


Change in pressure 


Pr 


Prandtl number 


q 


Rate of heat transfer 


Re 


Reynolds number 


AT j 
1 m 


Logarithmic mean temperature 


V 


Volume of exchanger 


w 


Mass flow rate 


e 


Exchanger heat transfer effectiveness 


P 


Mass density 


V 


Absolute viscosity 



Subscripts : 



A 


Material A 


B 


Material B 
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III. RESULTS 



In this study, the objective use to develop a series of 
computer programs to carry out. preliminary design on heat ex- 
changers required for thermal cycles. For this development, 
general properties of heat exchangers were examined and mathe- 
matical models were developed for various modes of heat trans- 
fer. 

The computer programs were developed and tested. All of 
the programs performed as desired except for the cost calcula- 
tions and the author feels that they should perform success- 
fully in studies of a broad range of inputs. The programs are 
explained and written in detail in the Appendix. 

Two of the programs were used in detail in an examination 
of the proposal made by Prof. Zener. See Aopendix G. The evan- 
oration and condensation heat exchanger programs carried out the 
desired preliminary design and showed the variety of designs 
that can be carried out with these programs. 

The Zener proposed cycle was examined to determine its 
feasibility. With the use of freon-21 as a working fluid, the 
evaporator must have a volume of two billion cubic feet for a 
100 Mega-watt power plant. The condenser will have a minimum 
size of over three-hundred million cubic feet of volume. More 
detailed results are included in Appendix G. 
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IV CONCLUSIONS 



The following conclusions can be drawn from the results of 
the study. 

1. The computer can be a valuable aid in heat exchanger 
design work . 

2. The cycle proposed by Prof. Zener is not feasible at 
the present time. 

3. The cost analysis done in this study is very inaccurate. 

4 . Although the programs can give a variety of geometrical 
arrangements there is no way to compare these arrange- 
ments to determine the best. 
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V RECOMMENDATIONS 



This study appears to the author to be weak in two areas . 
The evaluation of cost is very inaccurate. The problem of opti 
mization was not discussed in this study. A further study of 
these two areas would be beneficial. As these two areas are 
inter-related, they could be examined together. 

More examination of the cost of heat exchangers needs to 
be carried out. The relationshiDS for cost caluculation deri- 
ved in section 2.1 are based in I960 data. This in in itself 
is a problem, but also the size of heat exchangers designed is 
way beyond the range of validity of these relationships. To 
extrapolate to *4,000 feet a curve that is valid up to *'00 feet 
is unrealistic. More up to data values of cost should be de- 
termined. Also cost relationships that are valid for large 
heat exchangers should be derived. These relationships could 
easily be placed in the program. 

The programs as written contain no optimization operations 
The programs require a disigner to input various geometric 
arrangements and then compare the output to determine wich is 
best for his application. With up to data cost relationships, 
an optimization operation could be incoroorated into the pro- 
grams giving the designer the best heat exchanger within some 
specified limit. 
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Appendix A - COUNTER-FLOW CONDENSING HEAT EXCHANGER PROGRAM 

A.l Instructions for operation and Input 

This program is written in FORTRAN IV and can be run as is 
on the Interdata computer of the Joint Civil and Mechanical 
Engineering Computer Facility at M.I.T. With modifications in 
the control, input and output cards, this program can be run 
on any computer that is programed to compile FORTRAN IV.' 

For this program, the input data is broken into two groups. 
The first contains four cards, three of which contain properties 
of the fluids and the fourth is the number, M, of input sets to 
follow. The second group has 2M cards. There are two cards to 
a set. The first card serves as an identification card and the 
second has the geometric arrangement of the heat exchanger; Din, 
Dout , S, N, KMAT, and NB. With this input arrangement, M varia- 
tions of the geometry of the exchanger can be examined using the 
same fluid properties. 

The input data is punched on the cards as follows: 

(1) All values will be in units of lb , hr., 
feet, Btu , and °F. 

(2) All values will be fixed point numbers with 
decimal points except the value of M. 

(3) No commas are needed to separate the input 
values . 

(ty) The input values must be in the proper col- 
umns but need not be right or left justified. 

(5) All fixed point numbers can have up to 4 (four) 
digits after the decimal point. 

(6) The properties of the fluids, unless specified 
to be at a particular point, are average values 
over the length of the exchanger. 
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CARD 



COLUMNS 





1—13 


14—26 


27—39 


OJ 

in 

1 

1 

o 

-=r 


53—65 


66—78 


1 


TCI 


TC2 


Tw 


DELHW 


MKRW 


ROWL 


2 


ROW 


MUWL 


KWL 


KC 


PRC 


MUC 




1—12 


13—24 


25—36 


37—48 


49—60 




3 


ROC 


HFG 


CPC 


PCI 


PC2 




4 


1—4 












5 


Identification name 


for each 


geometric 


variation 






1—13 


14 — 26 


27—39 


40 — 52 


53—65 


66—78 


6 


DIN 


DOUT 


S 


N 


KMAT 


MB 


7 


Same as 


Card 5 










8 


Same as 


Card 6 











This repeats M times. 

The output will appear with the identification^at the top of 
a page and the following printed below it. 



MASS FLOW RATE WORKING FLUID 
MASS FLOW RATE COOLING FLUID 
NUMBER OF TUBES 
INNER DIAMETER OF TU3ES 
OUTER DIAMETER OF TUBES 
SPACING OF TUBES 
LENGTH OF TUBES 

DROP IN PRESSURE COOLING FLUID 
COMPARITIVE COSTS 



For each identification there will be a separate page. 
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A. 2 LIST OF VARIABLES 



AREA 

COST 

CPC 

CTB 

CW 

CV/V 

DELHV/ 

DIN 

DLPC 

DOUT 

DTLM 

GC 

HC 

HFG 

HWA 



KC 



KWL 



KM AT 



LENGTH 

M 



Heat transfer surface area in ft 2 

Cost of heat exchanger in dollars 

Specific heat at constant pressure of cooling 
fluid in Btu/lb m °F 

Cost of the materials for the tubes and shell 
in dollars 

Fabrication costs in dollars 
Fabrication cost per tube in dollars 
Change in enthalpy in Btu/lb^ 

Inner diameter in feet 

Frictional pressure drop in cooling fluid in 
lb/ft 2 

Outer diameter of tubes in feet 

Log mean temperature difference in °F 

Mass flow velocity in lb m /hr. ft. 2 

Heat transfer coefficient of cooling fluid in 
Btu/hr . ft. °F 

Latent heat of condensation in Btu/lb^ 

Average heat transfer coefficient of working 
fluid in Btu/hr. ft. 2 °F 

Thermal conductivity of cooling fluid in 
Btu/hr. ft. °F 

Thermal conductivity of working fluid as a 
liquid in Btu/hr. ft. °F 

Thermal conductivity of wall material in 
Btu/hr. ft. °F 

Length of heat exchanger in feet 
Number of groups of inputs data 
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A. 2 LIST OF VARIABLES (cont.) 



MFRC 

MFRW 

MUC 

MUWL 



N 

NB 

PCI 

PC2 

PRC 

Q 

REC 

ROC 

ROWL 

ROWV 

s 

TCI 

TC2 

TW 

TLF 

TMF 

U 



Mass flow rate of cooling fluid in lb m /hr. 

Mass flow rate of working fluid in lb m /hr. 

Viscosity of cooling fluid in lb m /hr. ft. 

Viscosity of working fluid as a liquid in 
Ib^/hr. ft. 

Number of tubes 

Number of tubes in a verticle row 
Inlet pressure of cooling fluid in lb/ft. 2 
Outlet pressure of cooling fluid in lb/ft 2 
Prandtl number of cooling fluid 
Rate of heat transfer in Btu/hr. 

Reynolds number of cooling fluid 

Density of cooling fluid in lb m /ft. 3 

Density of working fluid as a liquid in lb m /ft. 3 

Density of working fluid as a vapor in lbm/ft . 3 

Spacing of tubes 

Inlet temnerature of cooling fluid in °F 
Outlet temperature of cooling fluid in °F 
Temperature of working fluid in °F 
Tube length factor in dollars/ft. 

Tube material factor 

Overall heat transfer coefficient in BTU/hr.ft 2 °F. 
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NEAL N f i v f R w f RERC ,UB, RUM. ,KUL ,MUC , KN'AT , HW A , HFG ,HCt LENGTH, KC 
REAC(u, 1 CG 0 > TC 1 » TC 2 » Tfc » CEU«h»MFRW,PCM 
F t A L ( 3 , K ui ) R C I V , MUN L , Kk l. , KC , PR C » VU C 
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Appendix B - COUNTER-FLOW EVAPORATION HEAT EXCHANGER PROGRAM 

B.l Instructions for operation and input 

This program is designed to be run on the Interdata com- 
puter of the Joint Civil and Mechanical Engineering Computer 
Facility at M.I.T. The program is written in FORTRAN IV and 
can be run on any computer that is Drogramed to compile FOR- 
TRAN IV. If run on a computer other than the Interdata, some 
modification may be required to the control, input and output 
cards . 

The input data is broken into two sections. The first con- 
tains five cards, four contain fluid properties and the fifth 
is the number, M, of input sets to follow. The second section 
has 2M cards. The first card serves as an identification card 
and the second card has the geometric arrangement of the heat 
exchanger; DIN, DOUT, S, N, KMAT. With this input arrangement, 

M different variations of geometry can be examined using the 
same fluid properties. 

The input data is punched on cards with the following, 
format : 

(1) All values will be in units of lb m , hr, 
feet, Btu, and °F. 

(2) All values except M will be fixed point 
numbers and will have decimal points. 

M is a floating point number. 

(3) No commas are needed to separate the input 
values on a single card. 

00 The input values must be in the designated 
column but need not be right or left Justi- 
fied. 
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(5) AH fixed point numbers can have up 
to four digits after the decimal 
point. : 

(6) The fluid properties, unless speci- 
fied to be at a particular noint, 
are average values. 



CARD COLUMN 





1—14 


15—28 


29—42 


43—56 


57—72 


1 


DELHW 


TH1 


TH2 


TSAT 


ROWL 


2 


ROWV 


MUWL 


MUVW 


CPWL 


KWL 


3 


MFRW 


PHI 


PH 2 


PW1 


PW2 


4 


HFG 


CPH 


MUH 


ROH 


KH 




1—4 










5 


M 










6 


Identification name 


is columns 5 18 






1—14 










7 


DIN 


DOUT 


S 


N 


KM AT 


8 • 


Same as 


card 6 








9 


Same as 


card 7 








This is 


repaeated M times 








The 


output 


will appear 


with the 


identification 


name at 



top of the page and the follov.'ing printed below it. 

MASS FLOW RATE WORKING FLUID 
MASS PLOW RATE HEATING FLUID 
NUMBER OF TUBES 
INNER DIAMETER OF TUBES 
OUTER DIAMETER OF TUBES 
SPACING OF TUBES 



(cont. overleaf) 
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LENGTH OF TUBES 

DROP IN PRESSURE WORKING FLUID 
DROP IN PRESSURE HEATING FLUID 
COMPARITIVE COSTS 



For each identification there will be a separate page 



B. 2 LIST OF VARIABLES 



AREA 




Heat transfer surface area in ft? 


COST 




Cost of heat exchanger in dollars 


CPH 




Specific heat at constant pressure of heating 
fluid in Btu/lb m °F 


CPWL 




Specific heat at constant pressure of working 
fluid as a liquid in Btu/lb m °F 


CTB 




Cost of the material for the tubes and shell 
in dollars 


CV, T 




Cost of fabrication in dollars 


cwv 




Fabrication cost per tube in dollars 


DELHI-/ 




Change in ethalpy of v.’orking fluid in Btu/lb m 


DELPH 




Change in pressure of heating fluid, in lb/ft? 


DELPVJ 


(1) 


Frictional oressure drop of working fluid 
per segiment of exchanger in lb/ft? 


DELZ 


(1) 


Length of each segiment in feet. 


DELX 




Change in quality in each segiment 


DEQ 




Equivalent diameter in feet 


DIN 




Inner diameter of tubes in feet 


DOUT 




Outer diameter of tubes in feet 


FXTT 




Reynolds number factor 


GH 




Mass flow velocity of heating fluid in 
lb m /hr. ft? 


GW 




Mass flow velocity of working fluid in 
lb m /hr. ft. 3 


HH 




Heat transfer coefficient of heating fluid in 




. 


Btu/hr . ft 2 °F 


HFG 




Latent heat of evaporation in Btu/lb m 



-fi2- 



B. 2 LIST OP VARIABLES (cont.) 



HMAC 

HW 

IV 

KH 

KWL 

M 

MFRH 
MFRW 
MUH 
MUV7 L 



MUVW 

M 

PHI 

PH2 

PY/1 

PW2 

PRWL 

REWL 

ROH 

ROV/L 



Heat transfer coefficient for macroconventive 
effect in Btu/hr. ft 2 . °F 

Heat transfer coefficient for working fluid 
as a vapor is Btu/hr. ft 2 . °F 

Two phase frictional multiplier 

Thermal conductivity of heating fluid in 
Btu/hr. ft. °F 

Thermal conductivity of working fluid as a 
liquid in Btu/hr. ft. °F 

Number of groups of input data 

Pass flow rate of heating fluid in lbm/hr. 

Mass flow rate of working fluid in lb m /hr. 

Viscosity of heating fluid in lb m /hr. ft. 

Viscosity of working fluid as a liquid in 
lb m /hr. ft. 

Viscosity of working fluid as a vanor in 
lb m /hr. ft. 

Number of tubes in exchanger 

Inlet Dressure of heating fluid in lb/ft 2 . 

Outlet pressure of heating fluid in lb/ft 2 . 

Inlet pressure of working fluid in lb/ft 2 . 

Outlet pressure of working fluid in lb/ft 2 . 
Pradtl number of working fluid as a liquid 
Reynolds number of working fluid as a liquid 
Density of heating fluid in lb m /ft 3 . 

Density of working fluid as a liquid in lb m /ft 3 . 
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B.2 LIST OF VARIABLES (cont) 



ROWV 

REH 

S 

TH1 

TH2 

TLF 

TMP 

TSAT 

X 

XTT 

U 

Z 



Densitv of working fluid as a vapor in 
lb m /ft\ 

Reynolds number of heating fluid 
Spacing between tubes 

Inlet temperature of heating fluid in °F 

Outlet temperature of heating fluid in° F 

Tube length factor 

Tube material factor 

Saturation temperature in °F 

Quality 

Martinelli parameter 

Overall heat transfer coefficient in Btu/hr. ft 2 
Length of heat exchanger in ft. 
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Appendix C. COUNTER FLOW HEATER HEAT EXCHANGER PROGRAM 

C.l Instructions for Operation and Input 

This program is written in FORTRAN IV and can be run as is 
on the Interdata computer of the Joint Civil and Mechanical 
Engineering Computer Facility at M.I.T. The program can be run 
on any computer facility that is designed with the ability to com- 
pile FORTRAN IV. If run on a computer other than the Interdata, 
some modifications may be required to the control, input and out- 
put cards. 

The input data is broken into two sections. The first sec- 
tion contains eight input data cards. The first seven of these 
cards contain the fluid properties and the eight is the number, 

M, of input sets in the second section. The second section has 
2M cards. The first card in each of these M sets serves as an 
identification card. The second card has the data on the geometric 
arrangement of the heat exchanger; DIN, DOUT, and S. With this 
input arrangement, M different geometrical arrangements can be ex- 
amined using the same fluid properties. 

The input data is punched on cards as follows: 

(1) All values will be in units of lb m , hr., 
feet, Btu and °F. 

(2) All values except M will be fixed point 
numbers and will have a decimal point. M is 
a floating point number. 

( 3 ) No commas are needed to separate the input 
data on a single card. 
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(4) The input values must be in the desig- 
nated column but need not be right or left 
justified. 

(5) All fixed point numbers can have up to 4 
(four) digits after the decimal point. 

(6) The fluid properties, unless soecified 
to be at a particular point, are aver- 
age values. Any prooerty that the last 
figure of its variable name is a number 
is specified to be at that temperature. 



JARD COLUMN 





1—12 


13—2*1 


25—36 


37— *18 


*19—60 








1 


MFRW 


TH1 


TH2 


TV/1 


TV/ 2 








2 


PHI 


PH2 


PW1 


PW2 


MUK 








3 


MUM 


CPH 


CPW 


ROH 


ROW 








*1 


KH 


KW 


DCPW1 


DCPW2 


DCPV/3 










1—10 


11 — 20 


21—30 


31— 40 


*11—50 


51—60 


61—70 


71—80 


5 


DCPH1 


DCPH2 


DCPH3 


TV/ 3 


TH3 


DMUW1 


DMUW2 


DMUV/3 


6 


DMUH1 


DMUH2 


DMUH3 


DKW1 


DKV/2 


DKV/3 


DKH1 


DKH2 


7 


DKH3 


DROW1 


DROW2 


DROV/3 


DROH1 


DROH2 


DROH3 





1 — *1 
8 M 



5 18 

Identification name 
1—12 ' 13—2*1 25—36 

10 DIN DOUT S 

11 Same as 9 

12 Same as 10 
This repeats M times 



(cont. overleaf) 
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The output will appear with the identification name at the 



top of a page and the following printed below it. 



MASS FLOW RATE WORKING FLUID 
MASS FLOW RATE HEATING FLUID 
NUMBER OF TUBES 
INNER DIAMETER OF TUBES 
OUTER DIAMETER OF TUBES 
SPACING OF TUBES 
LENGTH OF TUBES 

DROP IN PRESSURE WORKING FLUID 
DROP IN PRESSURE KEATING FLUID 
COMPARITIVE COSTS 



For each identification there will be a separate page. 
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LIST OF VARIABLES 



AREA 

COST 

CPH 

CPU 

CTB 

CW 

CWV 

DIN 

DOUT 

DEQ 

DPH 

DP V/ 

DTLM 

D 1 

D 2 

D 3 

GH 

GW 

KH 

KW 

LENGTH 



Heat transfer area in ft 2 . 

Cost of heat exchanger in dollars 

Specific heat at constant Dressure of heating 
fluid in Btu/lb m °F 

Soecific heat at constant pressure of working 
fluid in Btu/lb m °F 

Cost of materials for shell and tubes in dollars 
Fabrication costs in dollars 
Fabrication costs per tube in dollars 
Inner diameter of tubes in ft. 

Outer diameter of tubes in ft. 

Equivalent diameter in ft. 

Frictional nressure dron in heating fluid in 
lb/ft 2 . 

Frictional pressure droD of working fluid in 
lb/ft 2 . 

Log mean temperature difference in °F 

Property at temperature 1 

Property at temnerature 2 

Property at temperature 3 

Fass flow velocity of heating fluid in lb m /hr. ft 2 

Mass flow velocity of working fluid in lb^/hr. ft 2 

Thermal conductivity of heating fluid in 
Btu/hr. ft. °F 

Thermal conductivity of working fluid in 
Btu/hr. ft. °F 

Length of heat exchanger in ft. 
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LIST OP VARIABLES (cont). 



MFRH 

MFRW 

MUH 

MUW 

M 

N 

PHI 

PH2 

PW1 

PW2 

PART 

PRH 

PRW 

REH 

REW 

ROH 

ROW 

S 

TH1 
TH2 
TH3 
TW1 
TV/ 2 
TV/ 3 



Mass flow rate of heating fluid in lb n /hr 

Mass flow rate of working fluid in lb m /hr. 

Viscosity of heating fluid in lb m /hr. ft. 
Viscosity of working fluid in lb n /hr. ft. 

Number of input data groups 
Number of tubes 

Inlet pressure of heating fluid in lb/ft 2 . 
Outlet pressure of heating fluid in lb/ft 2 . 
Inlet pressure of working fluid in lb/ft 2 . 
Outlet pressure of working fluid in lb/ft 2 . 
Temperature ratio 
Prandtl number of heating fluid 

Prandtl number of working fluid 

Reynolds number of heating fluid 
Reynolds number of ’working fluid 
Density of heating fluid in lb n /ft 3 . 

Density of working fluid in lb m /ft 3 . 

Tube spacing 

Inlet temperature of heating fluid in °P 
Outlet temperature of heating fluid in °F 
Intermediate temperature of heating fluid in °F 
Inlet temperature of working fluid in °F 
Outlet temperature of working fluid in °P 
Intermediate temperature of working fluid in °F 
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LIST OP VARIABLES 



( cont ) . 



TW1 Inlet temperature of working fluid in P 

TV. T 2 Outlet temperature of working fluid in P 

TW3 Intermediate temperature of working fluid in P 

TLP Tube length factor 

IMF Tube material factor 

XI X9 are grouped constants 

XXL XX12 are grouped constants 

X (J) Above properties for each segment of exchanger 
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Appendix D 



COUIITER-PLOW COOLER HEAT EXCHANGER PROGRAM 



D.l Instructions for Oneration and Inputs 

This program is v.'ritten in FORTRAN IV and can be run as is 
on the Interdata computer of the Joint Civil and Mechanical 'Engineering 
Computer Facility at M.I.T. The program can be run on any com- 
puter facility that is designed to compile FORTRAN IV. If run 
on a computer other than the Interdata, some modifications may 
be required to the control, input and output cards. 

The input deck is broken into two sections. The first con- 
tains eight input cards. The first seven cards contain the fluid 
properties and the eighth is the number, M, of input sets in the 
second section. The second section has 2M cards. The first card 
of each set serves as an identification card. The second card con- 
tains the data on the geometric arrangement of the heat exchanger; 

DIN, DOUT, and S. With this arrangement of the input data M differ- 
ent geometrical arrangement can be examined using the same fluid 
properties . 

The input data is punched on cards as follov.'s: 

(1) All values will be in units of lb^, hr., 
feet, Etu and °F 

(2) All values except M will be fixed point 
numbers and will have a decimal point. 

M is a floating point number 

(3) No commas are needed to separate the input 
data on a single card 

( ) The input values must be in the designated 
column but do not need to be right or left 
justified 
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(5) All fixed point numbers can have up to 

A (four) digits after the decimal point. 

(6) The fluid properties, unless snecified 
to be that at a particular point, are 
average values. Any property that the last 
figure of its variable name is a number is 
specified to be at that temperature. 



:ard column 





1—12 


13— 2A 


25—36 


37 — 48 


A9— 60 








1 


TCI 


TC2 


TW1 


TW2 


PW1 








2 


PW2 


PCI 


PC2 


MFRW 


MUW 








3 


CPW 


MUC 


CPC 


KW 


KC 








4 


ROW 


ROC 


DCPW1 


DCPVJ2 


DCPW3 










1—10 


11—20 


21—30 


31 — A0 


Al — 50 


51—60 


61 70 


71—80 


5 


DCPC1 


DCPC2 


DCPD3 


TW3 


TC3 


DMUW1 


DMUV/2 


D MU 17 3 


6 


DMUC1 


DHUC2 


DMUC3 


DKW1 


DKV/2 


DKW3 


DKC1 


DKC2 


7 


DKC3 


DR0W1 


DROV/2 


DR0V. T 3 


DR0C1 


DROU2 


DR0W3 





1 A 

8 M 



5 18 

Identification name 
1—12 13--2A 25—36 

DIN DOUT S 
•1 Same as 9 
Same as 10 

‘hese last two cards are repeated M times 



(cont. overleaf) 
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The output will appear with the identification name at 



the top of a page and the following printed below it. 



MASS FLOW RATE WORKING FLUID 
MASS FLOW RATE COOLING FLUID 
NUMBER OF TUBES 
INNER DIAMETER OF TUBES 
OUTER DIAMETER OF TUBES 
SPACING OF TUBES 
LENGTH OF TUBES 

DROP IN PRESSURE WORKING FLUID 
DROP IN PRESSURE COOLING FLUID 
COMPARITIVE COST 

For each identification name there will be a separate page. 
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LIST OF VARIABLES 



AREA 


Heat transfer area in ft 2 . 


CPC 


Specific heat at constant oressure for cooling 
fluid in Btu/lb m °F 


CPW 


Soecific heat at constant pressure of working 
fluid in Btu/lb m °F 


COST 


Cost of heat exchanger in dollars 


CTB 


Cost of materials for shell and tubes in dollars 


CW 


/ 

Fabrication costs in dollars 


CWV 


Fabrication costs per tube in dollars 


DIN 


Inner diameter of tubes in ft. 


DOUT 


Outer diameter of tubes in ft. 


DEQ 


Equivalent diameter in ft. 


DPC 


Frictional oressure drop of cooling fluid in 
lb/ft 2 . 


DPW 


Frictional pressure drop of working fluid in 
lb/ft 2 . 


DTLM 


Log mean temperature in °F 


D 1 


Property at temperature 1 


D 2 


Property at temperature 2 


D 3 


Property at temperature 3 


GC 


Mass flow velocity of cooling fluid in lb^/hr. ft 2 


GW 


Mass flow velocity of working fluid in lb m /hr. ft 2 


KC 


Thermal conductivity of cooling fluid in 
Btu/hr . ft. °F 


KW 


Thermal conductivity of working fluid in 
Btu/hr. ft. °F 


LENGTH 


Length of heat exchanger in ft. 
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LIST OF VARIABLES (cont) 



M 

MFRC 

MFRW 

MUC 

MUW 

N 

PCI 

PC2 

PW1 

PW2 

PRC 

PRVJ 

PART 

REC 

REW 

ROC 

ROW 

S 

TCI 

TC2 

TC3 

TLF 

TMF 



Number of inputs data groups 

Mass flow rate of cooling fluid in lb m /hr. 

Mass flow rate of working fluid in lb m /hr. 

Viscosity of cooling fluid in lb m /hr. ft. 
Viscosity of working fluid in lb m /hr. ft. 

Number of tubes in exchanger 
Inlet pressure of cooling fluid in lb/ft 2 . 
Outlet pressure of cooling fluid in lb/ft 2 . 
Inlet pressure of working fluid in lb/ft 2 . 
Outlet pressure of working fluid in lb/ft 2 . 
Prandtl number of cooling fluid 

Prandtl number of working fluid 

Temperature ratio 

Reynolds number of cooling fluid 
Reynolds number of working fluid 
Density of cooling fluid in lb m /f t 3 . 

Density of working fluid in lb m /ft 3 . 

Spacing of tubes 

Inlet tem.Derature of cooling fluid in °F 
Outlet temperature of cooling fluid in °F 
Intermediate temperature of cooling fluid in °F 
Tube length factor 
Tube material factor 



- 10 ')- 



LIST OF VARIABLES 


( cont ) . 


TW1 


Inlet temperature of working fluid in °F 


TV/ 2 


Outlet temperature of working fluid in °F 


TW3 


Intermediate temperature of working fluid in °F 


Xl X9 


are grouped constants 


XXI XX12 


are grouped constants 


X (J) 


Above properties for each segment of exchanger 
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Appendix E. CROSS-FLOW HEATER HEAT EXCHANGER PROGRAM 

E.l Instruction for Operation and Incuts 

Written in FORTRAN IV, this program is designed to be run 
on the Interdata cornDuter of the Joint Civil and Mechanical 
Engineering Computer Facility at M.I.T. With modifications in 
the control, input and output cards, this program can be run on 
any computer that is designed to compile FORTRAN IV. 

The first three cards of the input data are used for proper- 
ties of the fluids. Ten cards are then needed to input the ten 

j 

bt ten matrix F. The nest card contains a single number, M. 

This is the number of geometric arrangements that will be examined. 
The input cards that follow are in two card sets. The first of 
these contains an identification name, the second has the data on 
geometric arrangement, DIN, DOUT, S, N, and KMAT. With this in- 
put arrangement, M different geometric arrangements can be ex- 
amined with the same fluid properties. 

The input data is punched on the cards as follows: 

(1) All values will be in the units of 
lb m , hr., feet, Btu, and °F 

(2) All values exceot M will be fixed point 
numbers and will have a decimal Doint. 

M is a floating noint number. 

(3) No commas are needed to senarate the 
input data on a single card. 

(4) The input values must be in the desig- 
nated columns but do not need to be 
right or left Justified. 
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(5) All fixed point numbers can have up to 

4 (four) digits after the decimal point. 

( 6 ) The fluid properties are average values 
of these properties over the temperature 
range . 



CARD COLUMN 





1 — 12 


13 — 24 


25—36 


37 — 48 


£r 

'JD 

1 

CA 

o 


61—70 


1 


MFRW 


TWH 


TV/C 


TKH 


THC 


PHI 


2 


PII2 


PW1 


PW2 


CPW 


CPH 


MUV/ 


3 


MUH 


KW 


KH 


ROH 


ROW 





Cards 4 through 13 contain values of the matrix F 

Theses values are punched 10 to a card with six columns per value 

1 — 4 



14 


M 










5 18 








15 


Identification name 










1—12 13—24 


25—36 


37—48 


49—60 


16 


DIN DOUT 


S 


N 


KMAT 


17 


Same as 15 








18 


Same as 16 









These last two cards are repeated M times 

The output will appear with the identification name at the 
top of a page and the following printed below it: 

(cont. overleaf). 



-liii- 



MASS FLOW RATE WORKING FLUID 
MASS FLOW RATE HEATING FLUID 
NUMBER OF TUBES 
INNER DIAMETER OF TUBES 
OUTER DIAMETER OF TUBES 
SPACING OF TUBES 
LENGTH OF TUBES 

DROP IN PRESSURE WORKING FLUID 
DROP IN PRESSURE HEATING FLUID 
COMPARITIVE COSTS 



For each identification name there will be a separate page 



E . 2 LIST OF VARIABLES 



AREA 

COST 

CPH 

CPW 

CTB 

CW 

cwv 

DIN 

DOUT 

DLPH 

EC 

F 

GH 

GW 

HH 



HW 



KK 



KW 



LENGTH 

M 

MFRH 

*DTLH 



Heat transfer area in ft 2 . 

Cost of exchanger in dollars 

Specific heat at constant pressure of heating 
fluid in Btu/lb n °F 

Specific heat at constant oressure of working 

fluid in Btu/lb °F 

m 

Cost of materials for shell and tubes in dollars 
Fabrication costs in dollars 
Fabrication costs per tube in dollars 
Inner diameter of tubes in ft. 

Outer diameter of tubes in ft. 

Frictional pressure dron in lb/ft 2 . 

Heat exchanger effectiveness 
Mean temperature difference ratio 

Mass flow velocity of heating fluid in lb m /hr. ft 2 

Mass flow velocity of working fluid in lb m /hr. ft? 

Heat transfer coefficient of heating fluid in 
Btu/hr . ft 2 . OF 

Heat transfer coefficient of working fluid in 
Btu/hr. ft 2 . °F 

Thermal conductivity of heating fluid in 
Btu/hr. ft 2 , om 

Thermal conductivity of working fluid in 
Btu/hr. ft. °F 

Length of heat exchanger in ft. 

Number of input data grouns 

Mass flow rate of heating fluid in lb m /hr. 

Log mean ternerature difference in °F 
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LIST OP VARIABLES (E.2) (cont) 



MFRW 

MUH 

MUW 

N 

PHI 

PH2 

PW1 

PK2 

PRH 

PRW 

0 

REH 

REW 

P.OII 

ROW 

C 

o 

THC 

THH 

TVJC 

TWH 

TLF 

TMF 

U 

Z 



Mass flow rate of v/orking fluid in lb m /hr 
Viscosity of heating fluid in lb m /hr. ft. 

Viscosity of v/orking fluid in lb m /hr. ft. 

Number of tubes 

Inlet pressure of heating fluid in lb/ft 2 . 

Outlet pressure of heating fluid in lb/ft 2 . 

Inlet pressure of v/orking fluid in lb/ft 2 . 

Outlet pressure of v/orking fluid in lb/ft 2 . 

Prandtl number of heating fluid 
Prandtl number of v/orking fluid 
Heat transfer rate in Btu/hr. 

Reynolds number of heating fluid 
Reynolds number of v/orking fluid 
Density of heating fluid in lb m /ft 3 . 

Density of v/orking fluid in lb n /ft 3 . 

Tube spacing 

Cooler temperature of heating fluid in °F 

Hotter temnerature of heating fluid in °F 

Cooler temperature of v/orking fluid in °F 

Hotter temperature of v/orking fluid in °F 

Tube length factor 
Tube material factor 

Overall heat transfer coefficient in Btu/hr ft 2 °F 
Temperature coefficient 
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Appendix F. CROSS-PLOW COOLER HEAT EXCHANGER PROGRAM 



F.l Instructions for Operation and Inputs 

Written in FORTRAN IV, this program is designed to be run 
on the Interdata computer of the Joint Civil and Mechanical 
Engineering Computer Facility at M.I.T. With modifications in 
the control, input and output cards, this program can be run 
on any computer that is designed to compile FORTRAN IV. 

The first three cards of the input deck contain proper- 
ties for the fluids. The next ten cards are used to input the 
ten by ten matrix F. A single number, M, is placed on the four 
teenth card. This is the number of geometric arrangements that 
will be examined. The input cards that follow are In two card 
sets. The first of these two cards contains an identification 
name. The second contains the data on the geometrical arrange- 
ment of the exchanger, DOUT, DIN, S, N, and KMAT. With this in 
put arrangement, M different geometric arrangements can be ex- 
amined using the same fluid properties. 

The input data is punched on cards as follows: 

(1) All values will be in units of lb m , hr., 
feet, Btu, and °F 

(2) All values except M will have fixed ooint 
numbers and will have a decimal point. 

M is a floating point number 

( 3 ) No commas are needed to separate the 
Input data on a single card 

(*0 The input values must be In the designated 
column but need not be right or left justi- 
fied 



— 1 2 N — 



(5) All fixed point numbers can have ut> to 

4 (four) digits after the decimal point. 

(6) The fluid properties are average values 
averaged over the temperature range 

CARD COLUMN 





1—12 


13—24 


25—36 


37—48 


49—60 


61—72 


1 


MFRW 


TCH 


TCC 


TWH 


TWC 


PCI 


2 


PC2 


PW1 


PU2 


CPW 


CPC 


MUC 


3 


MUW 


KW 


KC 


ROW 


ROC 




Cards 


4 through 


13 contain values 


of the 


matrix F 




These 


values are 


■ punched 


10 to a 


card with 6 columns 


per value 




1—4 












14 


M 














5 


-18 










15 


Identification name 










1—12 


13—24 


25—36 


37—48 


49—60 




16 • 


DIM 


DOUT 


S 


N 


KMAT 




17 


Same as 


15 


- 








18 


Same as 


16 










These 


last two cards are 


repeated 


M times 







The output will appear with the identification name at the 
top of a page and the following printed below it: 



MASS FLOW RATE WORKING FLUID 
MASS FLOW RATE COOLING FLUID 
NUMBER OF TUBES 
INNER DIAMETER OF TUBES 



(cont. overleaf) 
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OUTER DIAMETER OP TUBES 
SPACING OP TUBES 
LENGTH OF TUBES 

DROP IN PRESSURE WORKING FLUID 
DROP IN PRESSURE COOLING FLUID 
COKPARITIVE COSTS 



For each identification name there will be a separate nage 



F . 2 LIST OF VARIABLES 



AREA 

COST 

CPC 

CPVI 

CTB 

CW 

CWV 

DIM 

DOUT 

DLPC 

DTLM 

EC 

F 

GC 

GW 

HC 

HW 



KC 



KVJ 

LENGTH 



Heat transfer surface area in ft 2 . 

Cost of exchanger in dollars 

Specific heat at constant pressure of cooling 
fluid in Btu/lb °F 

Specific heat at constant pressure of working 
fluid in Btu/lb m op 

Cost of materials for shell and tubes in dollars 

Fabrication costs in dollars 

Fabrication costs per tube in dollars 

Inner diameter of tubes in ft. 

Outer diameter of tubes in ft. 

Frictional pressure drop of cooling fluid in 
lb/ft 2 . 

Log mean temperature difference in °F 

Heat exchanger effectiveness 

Mean temperature difference factor 

Mass flow velocity of cooling fluid in lb m /hr. ft 2 

Mass flow velocity of working fluid in lb^/hr. ft 2 

Heat transfer coefficient of cooling fluid in 
Btu/hr . ft 2 . °F 

Heat transfer coefficient of working fluid In 
Btu/hr. ft 2 . °F 

Thermal conductivity of cooling fluid in 
Btu/hr. ft. °F 

Thermal conductivity of working fluid in 
Btu/hr. ft. °F 

Length of heat exchanger in ft. 
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F. 2 LIST OF VARIABLES (cont). 



M 

MFRC 

MFRW 

MUC 

MUW 

N 

PCI 
PC2 
PV/1 
PW2 
PRC 
PR V/ 

Q 

REC 

REV/ 

ROC 

ROV/ 

S 

TCC 

TCH 

TV/C 

TV/H 

TLF 



Number of input data groups 
Mass flow rate of cooling fluid in lb /hr. 

Mass flow rate of working fluid in lb m /hr. 

Viscosity of cooling fluid in lb m /ft. 
Viscosity of working fluid in lb m /ft. 
Number of tubes 

Inlet pressure of cooling fluid in lb/ft 2 . 
Outlet pressure of cooling fluid in lb/ft 2 
Inlet pressure of working fluid in lb/ft 2 . 
Outlet pressure of working fluid in lb/ft 2 
Prandtl number of cooling fluid 

Prandtl number of working fluid 

Heat transfer rate in Btu/hr. 

Reynolds number of cooling fluid 
Reynolds number of working fluid 
Density of cooling fluid in lb m /ft 3 . 
Density of working fluid in lb m /ft 3 . 

Tube spacing 



Cooler temperature of cooling fluid in F 
Hotter temperature of cooling fluid in °F 
Cooler temoerature of working fluid in °F 
Hotter temperature of working fluid in °F 
Tube length factor 
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P.2 

TMF 

U 

Z 



LIST OF VARIABLES (cont). 

Tube material factor 

Overall heat transfer coefficient in Btu/hr. ft 2 . °F 
Temperature coefficient 
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Appendix G 



EXAMINATION OP THE ZENER CYCLE 



G.l Examination 

To examine the usefulness of the computer progrmas de- 
veloped, the thermal cycle proposed by Prof. Clarence Zener(l) 
was examined. The program for designing condensers and the 
program for evaporator design were used. 

In his proposal, Zener suggests the use of ammonia as 
the working fluid. The condensor and the evaporator for this 
cycle were designed using ammonia and freon-21 as working fluids. 
Because of its high latent heat of vaporization about 5^0 Btu/lb m , 
ammonia requires a much larger heat exchanger for evaporation 
than does freon-21 with a latent heat of vaporization of 102 Btu/lb n . 
This difference in size can be seen in figures G-l, G-2 and G-3. 

The properties of Zener's prooosed cycle, as given In Ref. 1, 
were used as inputs to the computer programs. The allowable 
pressure drop over the length of the condenser and evaoorator 
was used as a constraint on the size of the exchangers. An 
inner diameter of \ inch was used as a starting diameter. If 
this diameter caused a frictional oressure drop greater than 
the allowed pressure drop, the diameter of the tubes was in- 
creased and the heat exchanger redesigned. 

In the condenser design, two values of the coefficient of 
heat transfer for condensation were used. The program was first 
run with the coefficient being calculated as discussed in section 



- 134 - 



2.3. The design was then revised using a coefficient of 
800 Btu/hr-ft 2 - °F. This is the value thatGregorig( 33) ob- 
tained using fluted tubes. This is further discussed in sec- 
tion 2.5. Zener in his proposal suggested the use of this tyoe 
of surface to increase the heat transfer coefficient. 

The two programs were run with the variation mentioned 
above and with different geometries. * The size of the heat ex- 
changers required for these different arrangements were then 
calculated . 

G.2 Results 

The condensor was sized using different numbers of tubes. 
The resultant required length of heat exchanger for each differ 
ent number of tubes are clotted in Figure G-l. As there is a 
significant improvement in reducing the size of the reauired 
condensor with the use of Gregorig surfaces, it is that size 
that is plotted. Plotted on Figure G-l are the results for 
both ammonia and freon-21. The minimum allowable diameter of 
the tubes for each given number of tubes is also plotted. 

Figure G-2 . 

The evaporator was sized in the same manner as the con- 
denser. Zener suggests the use of surface enhancement for the 
evaporator also. However, with the high mass flow rate surface 
enhancement appears to be of little advantage. The required 
length for the evaporator with the use of both ammonia and 
freon-21 is plotted in Figure G-3. 
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Minimum Diameter of Tubes in Feet 




Figure G-2 

Zener Cycle Condenser Humber of Tubes vers 
Tube Diameter 



- 137 - 



Length in Feet 



1500 



1250 ” 



100 $ 



7 50- 



5 00 



2 50 -—- 



0 . 



Ammonia 



25000 

Zener 



50000 
figure G-5 



75000 100000 

Humber of Tubes 



Cycle Evaporator Humber of Tubes vers 
Length 



- 138 - 



G. 3 Conclusions of the Zener Proposal 

The use of ammonia as a working fluid in this cycle apoears 
to be totally unrealistic. The minimum size requirement for 
the evaporator would amount to over thirty billion cubic feet 
with a length of over four thousand feet. The condenser using 
fluted tubes would be a great deal smaller, but still is over 
one billion cubic feet in total volume. 

Freon-21 is a much better working fluid for these heat ex- 
changers. The reduction in latent heat required to vaoorize 
the working fluid reduces the size requirement of the evaporator. 
Still the total volume required for the evaporator is two billion 
cubic feet. The condenser for a freon-21 system would be over 
three hundred million cubic feet in volume. 

The sizes of the above mentioned heat exchangers would 
be for production of 100 Mega-watts of electrical power. This 
is not an extremely large power plant by today's standards in 
terms of power output . 

From a thermodynamic point of view, this cycle is a feasible 
one. However, when the size of the required heat exchangers is 
taken into account, the cycle quickly appears to be unrealistic. 
V/ith a temperature difference of only two degrees Centigrade, in 
the evaporator, the mass flow rate is so great that the cycle 
is unusable. 
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